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We discuss the effects of so-called gluonic poles in twist-three hadronic matrix elements, as first 
considered by Qiu and Sterman, in the Drell-Yan process. These effects cannot be distinguished from 
those of time-reversal odd distribution functions, although time-reversal invariance is not broken by 
the presence of gluonic poles. Both gluonic poles and time-reversal odd distribution functions can 
lead to the same single spin asymmetries. We explicitly show the connection between gluonic poles 
and large distance gluon fields, identify the possible single spin asymmetries in the Drell-Yan process 
' and discuss the role of intrinsic transverse momentum of the partons. 
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■ I. INTRODUCTION 

On . 

In the usual description of the Drell-Yan process (DY) in terms of quark and antiquark distribution functions time- 
reversal symmetry implies the absence of single spin asymmetries at tree level, even including order l/Q corrections 
jjj. Additional time-reversal odd (T-odd) distribution functions are present when the incoming hadrons cannot be 
treated as plane- wave states. This may occur due to some factorization breaking mechanism |^|]. We will show that, 
i-C ■ even apart from such mechanisms, the contributions of T-odd distribution functions may effectively arise due to the 
presence of so-called gluonic poles attributed to asymptotic (large distance) gluon fields. The gluonic poles appearing 
in the twist-three hadronic matrix elements [^|~^) together with imaginary phases of hard subprocesses effectively give 
rise to the same single spin asymmetries. This is the origin of the single spin asymmetry of Ref. Q. Hence, the absence 
or presence of single spin asymmetries in DY can be viewed as a reflection of the absence or presence of gluonic poles. 
The "effective" T-odd functions coming from gluonic poles do not constitute a violation of time-reversal invariance 
itself. 

The outline of the article is as follows. We will first discuss how DY is described in terms of so-called correlation 
functions (section 2), which themselves are parametrized in terms of distribution functions (section 3). We focus 
especially on T-odd distribution functions, which show up in the imaginary part of the equations of motion (e.o.m.), 
which relate quark correlation functions with and without an additional gluon. In section 4 we will investigate the 
behavior of the quark-gluon correlation function in case it has a pole when the gluon has zero momentum. We will 
show that such poles will effectively contribute to the imaginary part of the e.o.m. and hence, to T-odd distribution 
functions. The large distance nature of gluonic poles is elaborated upon, in particular, the boundary conditions. In 
the final section we present the Qr-averaged DY cross-section with emphasis on the contributions of the effective 
T-odd distribution functions and the intrinsic transverse momentum. 



II. THE DRELL-YAN PROCESS IN TERMS OF CORRELATION FUNCTIONS 

We employ methods originating from Refs. ]s|-|l5|] in order to describe the soft (non-perturbative) parts of the 
scattering process in terms of correlation functions, which are (Fourier transforms of) hadronic matrix elements of 
non-local operators. We restrict ourselves to tree-level, but include l/Q power corrections. The asymmetries under 
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investigation are loosely referred to as 'twist-three' asymmetries, since they are suppressed by a factor of 1/Q, where 
the photon momentum q sets the scale Q, such that Q 2 — q 2 . We do not take Z bosons into account, since the 
asymmetries are likely to be negligible at or above the Z threshold. 

For the Drell-Yan process up to order 1/Q the quark correlation functions to consider are |g|- |ic| , |l^l 



$ ij (P 1 ,S 1 ;p) = J — I e*-*<P 1 ,S' 1 |V» J (0)^(«)|^,S' 1 ) ) (1) 

$% j (P 1 ,Sx;p 1 ,p2) = J ( ^)4 Q*y j Pl ' z e i<Pa - p&z '(Pi, SS^gA^z')^)^ Si). (2) 

We have included a color identity and g times t a from the hard into the soft parts <£> and <I>^, respectively. The 
inclusion of path-ordered exponentials, which arc needed in order to render the correlation functions gauge invariant, 
is implicit. 

The quark-quark correlation function can be expanded in a number of invariant amplitudes according to 

Dirac structure || . The available vectors are the momentum and spin vectors Pi , Si of the incoming hadron (spin- 
1/2), such that Pi -Si = 0, and the quark momentum p. In the case of a hard scattering process the momentum of the 
struck quark is predominantly along the direction of the hadron momentum, which itself is chosen to be predominantly 
along a light-like direction given by the vector n + . Another light-like direction ?i_ is chosen such that n + ■ n_ = 1; 
both vectors are dimensionless. The second hadron is chosen to be predominantly in the n_ direction, such that 
Pi ■ Pz = 0(Q 2 ). We make the following Sudakov decompositions: 



P 2 " = ^Z£ n '' + _^/ i (4) 



x 2 M? 



q ^^=n»+^=n»+q£, (5) 

for Qj, = —q\ <C Q 2 ■ We will often refer to the ± components of a momentum p, which are defined as p^ — p ■ n^. 
Furthermore, we decompose the parton momenta p,pi and the spin vector Si of hadron-one as 

_ xQ xi{p 2 +p 2 T ) 

:n + H -= — n_ +p T « xPi + pr, (6) 



x 



1V2 xQy/2 



pi= 7^ n + + n- +Pit ~ yP\ +P1T, (7) 

xiV2 yQy/2 

c - Al< 3 X1X1M1 Ai 

bi = -= n + -=- n_ + Sit ~ 7-ri + S 1T . (8) 

X1M1V2 QV2 Mi 

Also we note that up to order 1/Q only the transverse components of A a matter inside <&^. 

The Drell-Yan process consists of two soft parts and one of them is described by the above quark correlation 
functions, whereas the other is defined by the antiquark correlation functions, denoted by $ and <I>^. The correlation 
function $ depends on the second hadron momentum and polarization, P2 and S2, and the antiquark momentum k 
and is given by ||: 

1 ' j4„ -ik- 



<%(P 2 ,S 2 ;A0 = I d 4 z e- ife - z (F 2 ,^ 2 |^(z)^(0)|P 2 ,5 2 ). (9) 



The vectors in $ and are also decomposed in n±, 



,._ X Q 1 x 2 (k 2 + k 2 T ) ~-t,.u rm\ 

k — — n- H — — n + + kr ~ xP 2 -I- kx, (10) 



■f'2 



%/2 xQV2 



_ X2X2Q X2M2 A 2 

The function and the additional momentum ki are analogously defined as for the quark case. 
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At tree level four-momentum conservation fixes xP± = p + = q + = x\P^~ , i.e., x = X\ and similarly x = x%, 
and allows up to l/Q 2 corrections for integration over p~ and k + . However, the transverse momentum integrations 
cannot be separated, unless one integrates over the transverse momentum of the photon. In that case one arrives at 
correlation functions also integrated over their transverse momentum dependence, such that they only depend on the 
momentum fractions x,y and x,y. These partly integrated correlation functions <&(x), $(x), ^^{x, y) and ^^(x,y) 
are the quantities that are parametrized in terms of so-called distribution functions. For details see Ref. jjj. 

The five relevant diagrams lead to the following expression for the hadron tensor integrated over the transverse 
photon momentum (up to order l/Q): 



d 2 q T W"" = y |Tr ($(x)y>*$(x)-f u ) 

dy Tr ( n (y, xh»$(x)<y a j±= - ~ V . ^ ) + / dyTr ( *%(x, y)^^ - - f . -y a *(x)Y 



QV2 x-y + ie 'J J y ^ ' x - y + ie 

d * Tr {^^^^ih-'') -J dtTr (12) 

We will explain the above expression. The factor 1/3 arises from the color averaging in the qq annihilation. We have 
omitted flavor indices and summation; furthermore, there is a contribution from diagrams with reversed fermion flow, 
which is similar as the above expression but with fi <-> v and q — > — q replacements. In the expression the terms 
with fi± arise from the fermion propagators in the hard part neglecting contributions that will appear suppressed by 
powers of Q 2 , 

ix-i + rn _ ji+ x-y 



(pi - l) 2 — m2 + it QVZ x — y + ie 

$ — fti +m fi- x — y 

(q - ki) 2 - m 2 + ie Q^/2 x — y + ie' 



(14) 



where the approximate signs hold true only when the propagators are embedded in the diagrams. From these 
expressions one observes that the case x = y, i.e., the case of a zero-momentum gluon, corresponds to an on-shell 
quark propagator. 

Note that &(x), &(x), $>%(x, y) and ^^(x,y) are now integrals involving only one light-cone direction, for instance, 

^(x) = J ^e tXx (P,S\^(0)U^-)\P,S). (15) 

We observe that in the above expression one cannot simply replace ^^(x,y) by ^^(x, y), where 

*£«(s,tf)= J ^ge^e i "^){P ) ^|^(0)^(r ? n_)^(An_)|P, 1 9) (16) 

and iD a = id a + gA a . One must takes into account the difference proportional to J d 2 pTP^^{x,pT)- This difference 
is only zero, in case there are no transverse polarization vectors present. Similarly for the difference between and 

n- 

III. DISTRIBUTION FUNCTIONS 

For the correlation functions $ and $£ we need up to order l/Q the following parametrizations in terms of 
distribution functions |p^ , |l5| : 

$W = ^[/iWA +9i(x) Ai75 A +h 1 (x) l5 <$ 1T P 1 ] 



Mi 



e(x)l + g T {x)^fi 1T + /ii(a;)y75 [/£+, i-\ 



(17) 
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Mi 



Gd(i, y) i^Svrpfi + G D (i, y) 5f T7s A 



H D (x,y)X ll57 ^P 1 + E D (x,y)^P 1 



(18) 



where e T 



n+ a n-j3. We make a similar expansion for y) with the functions Gd, ■ ■ ■ replaced by Ga, 



while the rest stays the same. 

The parametrization of $(x) is consistent with requirements imposed on $ following from hermiticity, parity and 
time-reversal invariance, 



$t (Pi , Sup) = 70 $(Pi ,St\p) 70 [Hermiticity] 
$(Pi , 5i;p) = 70 $(Pi , -Si; p) 7 o [Parity] 
$*(Pi,5i;p) = 75 C$(Pi,5i;p) Ct 7s [Time reversal 



(19) 
(20) 
(21) 

where p = (p° 1 —p), etc. For the one-argument functions in Eq. ([Tt]) it follows from hermiticity that they are real. 
Note that for the validity of Eq. ( pl| ) it is essential that the incoming hadron is a plane wave state. For <& D and 
similarly for $J| hermiticity parity and time-reversal invariance yield the following relations: 



[$ D (P 1 ,S 1 ;p 1 ,p 2 ; Pi > Pi ) 7o [Hermiticity] 

$ D (Pi,Si;pi,p 2 ) = 7o$Da(Pi,-5 , i;pi,p 2 )7o [Parity] 
[$S(Pi ! 5i;pi,p 2 )]*=75C$ Da (Pi,5i;pi,p 2 )C t 75 [Time reversal] 

Hermiticity then gives for the two-argument functions in Eq. (|l^) the following constraints: 

Gd(x,v) = -G* D (y,x), 
Gd{x,v) = G* D (y,x), 
H D (x,y) = H* D {y,x), 
Ed{x,v) = -E* D (y,x). 



(22) 
(23) 
(24) 



(25) 
(26) 
(27) 
(28) 



Hence, the real and imaginary parts of these two-argument functions have definite symmetry properties under the 
interchange of the two arguments. If we would impose time-reversal invariance all four functions must be real and 
Gd and Hd are then symmetric and Gd and Ed are antisymmetric under interchange of the two arguments, such 
that at x — y only Gd and Hd survive. 

In the remainder of this section we do not impose time-reversal invariance and hence allow for imaginary parts of 
these functions. In addition, the following (T-odd) one-argument distribution functions then appear: 



<f(z)lT-odd =^ 



f T {x)e!^ 'S\TnlTv ~ e L (x)Xiij 5 + h(x)- fi. 



(29) 



Also we parametrize: 

J d 2 pTp^(x,p T ) = - ^ 



ifvr ( x ) it^SiTpPi - 9 [ rr( x ) S" t ^ 5 Pi 



a/3 



ht^(x)X ll5 ^P 1 +iht K1, (xh%P 1 



-L(l), 



The superscript (1) stands for the first fe^-momcnt of fey-dependent distribution functions /(a:, fc T ), 



d 2 k T 



(30) 



(31) 



This particular parametrization Eq. (30) is written in a form similar to Eq. (pLSj ), while using the fey-dependent 

functions of Ref. JlUfl , Note that fix (x) and h^ 1 ^ (x) are T-odd. 
We observe (since iD a — id a + gA a ): 
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dy [G D (x,y) + G D (y,x)] = / dy [G A (x, y) + G A (y, x)} - 2if^ 1} (x), 



dy 



dy 



IT 

Mi 



G A (x,y) + G A {y,x)\ + 2g\ 1 >{x), 
dy [H A (x,y)+H A (y,x)] - 2h^\x), 
dy[E A (x,y) + E A (y,x)}- 2ih^ 1] (x) , 



G D (x,y) + G D {y,x) 
dy [H D (x,y) + H D (y,x)] 
dy [E D {x,y) + E D (y,x)} 
while for the 'differences' no fe^-moments appear: 

J dy [G D (x,y)- G D (y, x)\ = J dy [G A (x, y) - G A (y, x)\ , 
etc. 

The two-argument functions and the one-argument functions are related by the classical e.o.m., which hold inside 
hadronic matrix elements jllj . Using the above parametrizations one has the following relations |l^,[l5| : 



(32) 
(33) 
(34) 
(35) 

(36) 



dy 
dy 



G D (x,y) - G D (y,x) + G D {x,y) + G D (y : x) = 2xg T (x) - 2—hi(x), 
G D (x,y) + G D (y,x) + G D (x,y) - G D (y,x) = 2ixf T (x), 



dy \ H D (x,y) +H D (y,x)^ = xh L (x) - —g^x), 



dy 



H D {x,y) - H D (y,x) 



-ixe L {x), 



dy \ E D (x,y) - E D (y,x)\ = xe(x) - T7h( x )> 



dy 



E D {x,y) + E D (y,x) = ixh{x) 



(37) 
(38) 
(39) 
(40) 
(41) 
(42) 



From this we see that the (T-odd) imaginary parts of the two-argument functions are related to the T-odd one- 
argument functions, as one expects. So if time-reversal invariance is imposed, the imaginary parts of the e.o.m. Eqs. 
(|38|), (40) and ( (42| ) become three trivial equalities. We like to point out that if one integrates Eqs. ( |37| ) and ( |38] ) over 
x, weighted with some test-function cr(x), one arrives at the sum rules discussed in 13 
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In order to observe the role of intrinsic transverse momentum, we will use some specific combinations of distribution 
functions, indicated by a tilde on the function. The tilde functions are the true interaction-dependent twist-three parts 
of subleading functions, which often contain twist-tw o p arts, (in analogy to #2) called Wandzura-Wilczek parts p8[ . 
They are defined such that in the analogues of Eqs. (|37|) - (j42|) for G A etc. only tilde functions appear, 



ReG A (x,y) + ReG A (x,y) = xg T {x) - Jj h i( x ) " Sit 0*0 = x 9t{x), 



1 it 
-C 1 )/ 



lmG A (x,y) +lmG A (x,y) = xf T (x) + f^\x) = xf T (x) 
dy 2 Re H A (x,y) 
dy 21m H A (x,y) 
dy 2 Re E A (x,y) 



xh L (x) - -r- r gi{x) + 2h,£ '{x) = xh L {x), 
ivi 

-xe L (x) ee -xe L (x), 



dy 2lmE A (x,y) 



xe(x) - -Tihix) = xe(x), 
xh(x) + 2h^ (1 \x) ee xh(x). 



(43) 
(44) 
(45) 
(46) 
(47) 
(48) 
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IV. GLUONIC POLES AND TIME-REVERSAL ODD BEHAVIOR 



We are interested in the behavior of the quark-gluon correlation function <£>^ in case x = y, when the gluon has 
zero-momentum. For this purpose, we define (a is a transverse index) 

$ a Fij {x,y) = |^g e ^e^"^(P,5|^(0)F+"(^_)^(An_)|P,5) (49) 

and F f " T {z) = jj [D p (z), D a (z)]. Defined as given above, the matrix element has the same hermiticity, but the 
opposite time-reversal behavior as and $^ and we will parametrize it identically with help of functions called 
Gp(x,y),Gp(x,y), Hp(x,y) and Ep(x,y), noting that time-reversal implies (in contrast to or $^) that Gf and 
Ep are symmetric and thus may survive at x = y. In the gauge A + — one has F +a = d + A^ and one finds by 
partial integration 

(x-y)*%(x,v) = -&F(x,v)- (50) 
If a specific Dirac projection of (f>^(a;,a;) is nonvanishing, then the corresponding projection of $>°h(x,x) has a 
pole, hence the name gluonic pole. An example is the function T(x,St) = ttTt $>p(x, x) erpaS T fi- / P + = 

2-KiMS^Gp{x,x) discussed by Qiu and Sterman in Ref. ||f||. 

In order to define Eq. (|5^) at the pole, one needs a prescription, which is related to the choice of boundary conditions 
on A a {rj — ±oo) inside matrix elements. Possible inversions of F +a — d + A[f are (only considering the dependence 
on the minus component): 

/oo 
dz- 0(z- -y-)F +a (z-) 
-CO 

/OO 
dz- 9(y- - z-)F +a {z-) 
-oo 

g gM^|M _i jy_ <z _ y _ )p ^ z _ y (5l) 

One can use the representations for the 9 and e functions, to obtain 

n& V) = Six - y) *5(oo)(») + x _~y + ie f) 

= 6(x - y) <^(_oo)(z) + x _~y_ i£ ®f(x, V) 

= S (x - y) ' - ' + P — V), (52) 

2 x - y 



S(x-y)^ A{±oo)lJ (x)^ I ^p-e i ^e^»-*HP,S\i> j (0)gA%(r, = ±oc)M*n-)\P,S). (53) 



where 

S(x - v) , - .(x) = I 

2tt 2tt 

So Eq. (|52| ) shows the importance of boundary conditions in the inversion of Eq. (|5Q|) , if matrix elements containing 
A™ (77 = ±00) do not vanish. When such matrix elements vanish (implicitly assumed in |Q) the pole prescription does 
not matter. Also one obtains 

2n<S> F (x 1 x) = [$2 (oo) (z) - n(-co)(x)} , (54) 

which shows the relation between the zero-momentum quark-gluon correlation function and the boundary conditions. 
The behavior of ^2(±oo) ( x ) un< 4er time-reversal is: 

®T(±oc)(x) =7 5 C^ oo)oi (x)C^ 5 . (55) 
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This relation implies that time-reversal invariance only allows for symmetric or antisymmetric boundary conditions. 

To study the effect of gluonic poles we will consider the (nonvanishing) antisymmetric boundary condition^], 
®1(oc)( x ) ~ ~^A(-oo)( x )> wn i cn implies n & F (x, x) = ^a(oc)( x )- ^ n the diagrammatic calculation resulting in Eq. 
( |l2| ) one always encounters the pole of the matrix element (in this case in the principal value prescription) multiplied 
with the propagator in the hard subprocess (having a causal prescription), 



a cff 



(y,x) = 



x — y + ie 



§ F {y,x)=$'X(y,x)-Tz8{x-y) ® F (y,x). 



(56) 



The time-reversal constraint applied to ^(x, y) implies the analogue of Eq. (pjj), while § F (x,y) has the opposite 
behavior under time-reversal compared to ^ A (x,y). Thus for $^ eff (x, y) one does not have definite behavior under 
T-reversal symmetry. Specifically, the allowed T-even functions of & F (x, x), Gf(x, x) and Ep{x, x), can be identified 
with T-odd functions in the effective correlation function (f>^ cff (x, y). This implies that G° A s (x,y) and E A s (x,y) will 
have an imaginary part and this gives rise to two "effective" time- reversal-odd distribution functions f T s (x) and 
h eS (x) via the (imaginary part of the) e.o.m. Since by identification 



Gp(x,x) = J dy!mGf(y,x), 



iir Ep{x, x) 



dylmEf{y,x), 



it follows from Eqs. (|4J) and (|48|) that 

xf T S {x) = inGpfa, x) - 

xh eS (x) = 2mE F (x,x) 



1 



2MS T 
—iir 



IMP 



T(x,S T ), 
: Tr[$p(x,x)j Ta ', 



(57) 
(58) 

(59) 
(60) 



The function receives no gluonic pole contribution, since time-reversal symmetry requires Hp(x,x) = 0. 

Of course, the mechanism for generating finite projections of <S> F {x,x) remains unknown. We just can conclude 
that if there is indeed a non-zero gluonic pole (in the case of non-zero antisymmetric boundary conditions), then at 
twist-three there are two non-zero "effective" T-odd distribution functions, namely fx and h. The first one generates 
the single spin twist-three asymmetry found by Hammon et al. in their notation it is proportional to T(x, x). The 
second one leads to a new asymmetry (see next section). Summarizing, we find for the parametrization of &a(oo)( x )' 



)(*)=" 



ixM 



ff(x) i efS Tp P+h cfi (x)^P 



(61) 



which is constrained by time-reversal symmetry but behaves exactly opposite to for instance $g(a;), hence in their 
parametrizations the meaning of time-reversal even or odd functions are opposite also. 

The case of nonvanishing symmetric boundary conditions is less interesting, since <f> F (x,x) = 0, but it is allowed. 
The delta-function singularity in this case will contribute to the functions Ga(x, x) and Ha(x, x) and hence, to T-even 
tilde functions. This would only affect the magnitude of (time-reversal even) double spin asymmetries. 

The antisymmetric nonvanishing boundary condition for ^a(±oo)( x ^ m ight arise from a linear A-field, giving a 
constant field strength (cf. e.g. |2(],|2lJ] ) . One might also think of an instanton background field. In both cases one 
should interpret infinity to mean 'outside the proton radius'. Also, the constant field strength should be understood 
as an average value of the gluonic chromomagnetic field, which is non-zero due to a correlation with the direction of 
the proton spin. The large distance origin of the asymmetries arising from such a gluonic pole is apparent. 

We like to point out that so-called fermionic poles play a role in off-forward scattering, such as prompt photon 
production ]2g|j-^|, but not in DY to this order. 

The fragmentation function that is the analogue of the distribution function fx (called Dx), shows up in a single 
spin asymmetry in hadron production in e + e _ annihilation ]23| , |2"4| ] , allowed because final state interactions lead to 
T-odd fragmentation functions. In Ref. p5| both gluonic poles and final state interactions are considered, but without 
taking into account boundary terms in the matrix elements. This result is in fact an example of the effective relation 
we have shown (see also [0). 



1 The consistency of antisymmetric boundary conditions with Maxwell's equations has already been shown in Jl8| |. 
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V. THE DRELL-YAN CROSS-SECTION IN TERMS OF DISTRIBUTION FUNCTIONS 



We will now discuss the Drell-Yan cross-section in case one integrates over transverse photon momentum. So one 
uses the above parametrizations of the correlation functions in the expression for the integrated hadron tensor as given 
in Eq. ([l2]), which after contraction with the lepton tensor yields the cross-section. The parametrizations in terms 
of distribution functions are defined with help of the vectors n+, n_ and several transverse vectors. However, angles 
we are going to discuss with respect to another set of vectors. Depending on the choice of this set, we find different 
combinations of functions with and without a tilde. Needless to say that the cross-section itself is an observable and 
does not depend on the choice of vectors, even though its appearance changes. 

We choose the following sets of normalized vectors: 



t = q/Q, 

x = qr/Qr = {q-X\P\- x 2 P 2 )/Qt 

characterized by a parameter c and where Pi = Pi — q/(2xi), such that: 

fj 1 
+ V2 

nl = —= 
V2 



*" + 



2c ~q x . 



- - 2(1 - c) 



Q 



So the parameter c basically distributes the transverse momentum between Pi and Pi in different ways (Fig 
c = (c — 1), then Pi (P 2 ) has no transverse component. The symmetric case c = 1/2 is the one used in Ref. [2 



(62) 
(63) 
(64) 

(65) 

(66) 
]). If 




lepton plane (cm) 

FIG. 1. Kinematics of the Drell-Yan process in the lepton center of mass frame, for a particular value of c. 



In this way we arrive at the following expression for the Drell-Yan cross-section in case of unpolarized leptons: 
da(hih 2 ^UX) a 2 ^ 2 J . - , 

a. a \_ 

+ B(y) \S 1T \ \S 2T \ cos(0 Sl +^s 2 ) (hihi) 

+ C(y) \S 1T \ sin(0 Sl ) (^xi ((1 - c)f T + c/r) 7 X + ^ x 2 h x (ch+(l-c)t 



C(y) \S 2T \ sin(<fe 2 ) x 2 f x (cJ T + (1 - c)%) ((1 - c)h + d, ) /, , 



+ C(y) X 2 \S 1T \ cos((f)s x ) x i (C 1 ~ c )9t + cg T )g 1 + '~^ L X2 h i [chh + (1 - c)hi 

-C(y) Ai \S 2T \ cos(05 2 ) (^^x 2 gi (cg T + (1 - c)g T ) + x i {i 1 - c ) h L + chtjhi 



(67) 



+ z^z v ^ l v , the perpendicular part of the 



where dQ = 2dy d(f> 1 and <p l gives the orientation of = (g^" — i^i"^ 
lepton momentum I, and y = l~ /q~. In this result we encounter the following functions of y: 

A(y) = (l-2y + 2y 2 ) /2, 

B(y) = y(l-y), 

C(y) = (1 - 2y) y/y{l-y). 



(68) 
(69) 
(70) 



Furthermore, fifi — f\{%\)f\{%2)i etc. and where a is the flavor index. 

For c = 1/2 we find agreement with the results of JTJ for the cross-section without T-odd distribution functions. 
Hence, we confirm the deviation of that result from the one found in [p"5| . 

We observe single-transverse-spin asymmetries with two possible angular dependences, namely sin(0s 1 ) and 
sin(</>s 2 ). Each of them comes with two products of functions, in particular an unpolarized one (fx or h) times a 
polarized one (/t or h\). There is no choice of c to eliminate the tilde functions from this expression, nor to only 
retain tilde functions. This shows the non-trivial role of intrinsic transverse momentum of the partons and one cannot 
discard it. This means that unlike in the case of DIS, one cannot take only $(x) and ^^(x,y) as a basis Jl2[ . 

If we assume that the presence of T-odd distribution functions is only effective, arising due to gluonic poles, and 



that $ 



A(oo) 



$« (oo)) then ff = ff and h 



cff 



This implies the following single spin asymmetry (hadron-two 



unpolarized), given in the lepton center of mass frame: 



A T = 



2sin(26>) sin( 



1 



'IT 



Q 



2M 1 x x f T (x 1 )f^(x 2 ) + 2M 2 h a l {x l )x 2 h a (x 2 ) 



elft{xi)tf{x2), (71) 



where we used that y — (1 + cos 9)/ 2 and 9 is the angle of hadron-two with respect to the momentum of the outgoing 
leptons. The first term in the asymmetry (proportional to fx) is the one discussed in JjJ (in their notation it is 
proportional to T(x,x)q(y)), which will also be present in DIS (fi(x 2 ) = 5(1 — x 2 )). The second term is the other, 
new single spin asymmetry arising in DY from a gluonic pole. It is not proportional to T(x,St), but to another 
projection of in the point x = y, cf. Eq. (60). 



VI. CONCLUSIONS 



We have shown how the effects of so-called gluonic poles in twist-three hadronic matrix elements, which were 
first discussed by Qiu and Sterman cannot be distinguished from that of T-odd distribution functions. We 

investigated this for the Drell-Yan process, which is expressed in terms of products of distribution functions. Even in 
the absence of T-odd distribution functions, imaginary phases arising from hard subprocesses together with gluonic 
poles give rise to effective T-odd distribution functions. This leads to single spin asymmetries for the Drell-Yan 
process, such as the one found recently by Hammon et al. Q. These asymmetries therefore can have a different origin 
than the analogous asymmetries in inclusive hadron production in e + e~ annihilation | 
to final state interactions, which are expected to be present always in contrast to initial state interactions. We have 
moreover shown that the presence of gluonic poles is in accordance with time-reversal invariance and requires a large 
distance gluonic field with antisymmetric boundary conditions. Our analysis shows also the role of intrinsic transverse 
momentum of the partons for the DY cross-section at subleading order. 



23 24], which can also arise due 
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